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Abstract. 

Present models describing the interaction of quantum Maxwell and gravitational 
fields predict a breakdown of Lorentz invariance and a non standard dispersion relation 
in the semiclassical approximation. Comparison with observational data however, 
does not support their predictions. In this work we introduce a different set of ab 
initio assumptions in the canonical approach, namely that the homogeneous Maxwell 
equations are valid in the semiclassical approximation, and find that the resulting field 
equations are Lorentz invariant in the semiclassical limit. 

We also include a phenomenological analysis of possible effects on the propagation 
of light, and their dependence on energy, in a cosmological context. 
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1. Introduction 

Present models describing the interaction of quantum Maxwell and gravitational fields 
predict a breakdown of Lorentz invariance. In the semiclassical approximation, the 
common feature in these models is a non standard dispersion relation which shows 
that the spacetime behaves as a media with a frequency dependent index of refraction 
121 HH H| - Another consequence is the selection of a preferred reference frame, namely, 
the one at rest with respect to the media. In geometrical terms this can be realized by 
the introduction of a preferred timelike vector field t a that serves as a universal time. 

It is instructive to analyze in pure phenomenologial terms the effect of a 
semiclassical gravity state on the propagation of electromagnetic radiation. The effective 
interaction is realized by adding to the standard Maxwell equations extra terms with 
minimal coupling between the Maxwell field and the above mentioned timelike vector 
field. Although the interaction contains parity breaking and parity preserving terms, it 
is perhaps surprising that the only possible mesurable effect arising from our approach 
comes from the parity breaking interaction and it is this effect, precisely, what is 
predicted by the leading models on the canonical approach [21 01 Ej. The main prediction 
of these models is a dispersion relation that depends on the helicity and the energy of 
the radiation field. Our analysis also seems to rule out wave propagation that is parity 
invariant, and in particular recent results by John Ellis, et. al j^]. 

Note that if the dispersion relation depends on the helicity then a linearly polarized 
wave packet with a continuum spectra will rotate its polarization direction as it travels 
through space. Since the rotation depends on the energy, the wave will become totally 
depolarized after traveling a suitable optical path. For light coming from a cosmological 
source, the observation of linear polarization can be used to set a severe bound on the 
phenomenological coupling constant of the Gambini-Pullin model [B] and, as we will 
show in this work, essentially rule out the Sahlmann- Thiemann construction. 

At this point one is tempted to ask whether or not Lorentz invariance is necessarily 
broken by quantum gravity. Since observational evidence seems to support Lorentz 
invariance, we reexamine the canonical models looking for assumptions introduced in 
their construction that, although natural at first sight, might not be true in the final 
theory. Our hope is that by suitably changing them we may preserve Lorentz invariance. 

In particular, in what follows we replace their ab-initio assumption of the electric 
field and potential being conjugate variables by a different anzatz, namely, that the 
quantum Maxwell field is the curvature of the Maxwell connection. This natural 
assumption has the expected consequence that the resulting field equations are Lorentz 
invariant. Although the calculation is done at a linear level, it is very easy to extend the 
results to any order in the perturbation parameter. We thus conjecture that Lorentz 
invariance will still hold in the quantum interaction between gravitons and photons. 

The work is organized as follows. In Section 2 we present a phenomenological 
approach to the propagation equations showing that if Lorentz invariance is broken, 
at most there will be a gravity induced rotation of the polarization vector. We also 
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show that the available observational data essentialy rules out the Sahlmann- Thiemann 
model. 

In Section 3 we assume that the Maxwell field is the curvature of a connection, 
namely, F = dA. Since this statement is independent of the existence of a quantum 
gravitational field we conjecture that it is valid for the full theory. Using the energy 
density of the electromagnetic field and this new assumption, we obtain different 
conjugate variables and thus different field equations for the Maxwell field from those 
obtained in previous canonical models. 

Finally, in the conclusions we summarize the main results obtained in this work. Is 
Lorentz invariance broken when light propagates on a quantum space-time? We address 
this question an add our own bias on the subject. 



2. Phenomenological equations for the propagation of light 

The construction of equations for the propagation of light under the influence of effects 
arising from the quantum nature of gravity requires explicit assumptions about the form 
that we expect the low energy regime of these effects will take. In the particular case of 
cosmological applications, we would expect that the quantum expectation values should 
be in consonance with the properties of the classical metric g a b, defining the classical 
geometry. For the standard models, this is characterized by a timelike vector field t a , 
whose integral lines are the world lines of comoving observers with 4-velocity u a . We 
take this as implying that the expectation values defining the (local) low energy limit 
should be tensor functions only of u a , g a b and a scalar function of t, where t is an affine 
parameter ("time") for the integral lines of t a . 

In view of some recent theoretical results pQ we take the Maxwell tensor F ab as the 
fundamental physical quantity describing the electromagnetic field. Its components in 
a local Lorentz frame are F 0i = and F^ = e^B^^ where Ej_ and are, respectively, 
the components of the electric and magnetic field vectors, and where is the Levi- 
Civita symbol. We, therefore, do not assume the existence of a vector potential A a . 
In the absence of quantum gravity effects, the propagation of light is governed by the 
equations for the electromagnetic field F ab , which may be written as 



V a F\ b = 0, (1) 
V a F ab = 4vrJ b , (2) 

where F* ab is the dual of F ab , J a is the electric current, and we recall that the vanishing 
of the right hand side of corresponds to the absence of magnetic type currents and 
monopoles. 

The modification of these equations that we envision is the addition of terms 
on the right hand sides of (JTJ) and (j2j), corresponding to the presence of different 
"effective currents" possibly resulting from quantum gravity effects and to which the 
electromagnetic field gets coupled. We shall restrict to terms that are linear in the 
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electromagnetic field, assuming that non linear effects are of higher order and may be 
neglected in this approximation. Our phenomenological equations take the form, 



V a F\ b = ^t a F ab + ij 2 t a F* ab 

+ ^t a t c W c F ab + ^ A t a t c W c F\ b 
+ ^ 5 t a t c t d V c [V d F ab ] + At a t c t d V c [V d F\ b \ 
+ *P 7 t a g cd V c [V d F ab ] + iP s t a g cd V c [V d F* ab }, (3) 
V a F ab = 4nJ b + X it a F ab + X 2t a F* ab 

+ X 3t a t c V c F ab + X it a t c V c F* ab 

+ X5 t a t c t d V c [V d F ab ] + X 6t a t c t d V c [V d F* ab ] 

+ Xit a g cd V c [V d F ab ] + X st a g cd V c [V d F* ab ]. (4) 

In these equations we are assuming local couplings that lead to expressions in F ab and 
its derivatives, and we have included only terms up to second derivatives. Notice that, 
from our assumption of a cosmological metric we have, 



t a V a t b = 0, 

Similarly, from simple physical arguments, we expect the coefficients x% an d i>i to 
depend only on the scale parameter ("radius of the Universe") of the metric, in such a 
way that they vanish when the Planck length £p is taken equal to zero. Note that even 
if some or all of the coefficients Xi an d ipi are non zero, on observational grounds they 
must be small. Thus, we may consider the effect of each term separately, as any cross 
effects would be of higher order. This analysis is carried out in the next subsection. 



2.1. Local Lorentz frame analysis 

The equations are more easily analyzed in a local Lorentz frame, adapted to the 
symmetry of the metric. Namely, if the local coordinates are (x,y,z,t), we have 
t a = (0,0,0, 1). As a first approximation we also neglect curvature effects, and equate 
covariant to ordinary partial derivatives. Moreover, we assume J a = 0, corresponding 
to "vacuum" propagation, and take all ipi, and Xi as constants. 

Using the standard forms for F ab and F* ab , we immediately find 



V-E = 0, 
V-B = 0, 



and 



(6) 
(7) 



VxE + d t B = X iE + ipiB + X 2d t E + ^ 2 d t B 
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+ X 3%E + ^d 2 t B + X4 V 2 £ + ij A V 2 B, 
VxB-d t E = X5 B + ^ b E + xedtB + ^d t E 

+ xrdtB + ij 7 d 2 E + Xs V 2 5 + ^V 2 E. 



(8) 
(9) 



We remark that in (JHJ) and © we have regrouped some terms from (JHJ) and (J1J), 
and renamed some constants. Notice that the factors of Xi violate parity conservation. 
It can also be seen from (JHJ) that the constants ipi, . . . , ^4, and Xi> ■ ■ ■ > X4 should vanish 
if one assumes absence of magnetic currents. In this latter case F a b admits a vector 
potential A a . This situation will be analyzed in the next Section. 



2.2. Plane waves 

We look now for plane wave solutions traveling along the z-axis. In this case, on account 
of ©, we should have, 



— * 

E = Re {(E x e x + E y e y ) exp(i(ut — kz))} , 
B = Re {(B x e x + B y e y ) exp(i(ut — kz))} . 



(10) 



Replacing (JTUJ) in (JHJ) and (JHJ) we find the general form for the dispersion relation 
k = k{u),Xi,ipi)- In view of the smallness of x% an d ipi, this may be expanded up to 
linear order in this coefficients, but it will be convenient to keep higher order terms. If 
we consider separately terms by their order of derivatives, and their parity conserving 
or violating character we find, in the parity conserving cases: 



k = \Ju 2 + 2iu(tpi + ^5) - 4^i^ 5 
~ uj + i(ipi +V5), 

k = Uy/{l- 2(^2 + + 4 ^ 6 ) 
~ UJ[1 - {lp 2 + A% 

k = uj\J[1 - 2iu(4> s + ip-j) - Au 2 i/j 3 i/j 7 ] 
~o;[l - iuj(ip 3 + ip 7 )}, 

k = uj yjl — k 2 (ip 4 - ip 8 ) 2 + ik(ip4 + ^8 
~ u[l - iuj(ip 4 + Vs)], 
and for the parity violating cases: 



k = sJuj 2 ± 2k(xi - Xs) + 4 XiX5 
~w±(xi-X5), 



\/l - (X2 + Xe) 2 ± *(%2 - Xe) 
w[l±i(x2 -Xe)], 



k = u) 



^/1+UJ 2 (X3 + X7) 2 ± ^(X3 - Xr) 



(11) 

(12) 
(13) 

(14) 



(15) 
(16) 
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- w[l±w(x3 - Xr)], 

k = u [l ± 2/c(x4 - Xs) - 4/c 2 X4Xs 

- w[l±w(x4 - Xs)]- 



1/2 



(17) 



(18) 



If we use these results to obtain the corresponding amplitudes, we find that in the 
parity conserving cases there is a single mode of propagation, and the polarization of 
a plane polarized wave is preserved in time, while in the parity violating cases there 
are two modes, one corresponding to right, and the other the left circular polarization. 
Possible observational effects for these cases have been discussed in jS]- 

Considering now particular cases, we see that (fTTjl corresponds to a frequency 
independent amplification or attenuation (depending on the sign of ipi + 0s) of the wave 
amplitude with time, with no effects on the polarization. The effect would be absent 
for ipx = —ips, but this would still leave a term quadratic in ipi, that would behave as 
a tachyonic mass added to the photon. The important question here is the order of 
magnitude of these couplings. In the absence of a theory we may only conjecture that 
we would expect these to be of the order of the corresponding dimensional quantities 
that characterize the model. In our case these are the Planck length £p, and time 
tp and possibly, for instance, the radius (or scale) a(T), and age T (or the Hubble 
constant H = a /a) of the Universe. Then, since ip\ and ^5 should have dimension 
[length] -1 , one would expect these quantities to be of the order of tpH 2 , which, after 
multiplication by the time to travel through cosmological distances, is too small to have 
any phenomenological relevance. 

In the case (II 2j) we have a frequency independent change in the velocity of the 
waves, which would be larger or smaller than that of light, depending on the sign of 
^2 +'06- Since ip2 and ipe are dimensionless, an estimate for them could be tpH, which 
is again too small for observable consequences. 

The case (J14jl is conceptually similar to (jllj) . leading again to amplification or 
attenuation of the waves. However, in this case ^4 and ^8 have dimension of [length]. 
Then, if we take them to be of the order of £p, the relevant quantity would be of 
order £pL/X 2 , where A is wavelength of the wave and L of the order of the distance 
to cosmological sources. Perhaps surprisingly, if we assume A = 10 _5 cm, (for visible 
light), and L = 10 9 ly, we find already £pL/X 2 ~ 10 3 . This would imply that visible 
light would not reach us but, in fact, not only visible light but also 7-rays are observed 
from cosmological sources. This implies that either 04 and ^8 are much smaller than 
this scale or 04 ~ —ips- If we take 04 = —"08! we find k ~ to — 2ip^u! 3 . This implies 
a group velocity v g — 1 — 6?/>f u; 2 which, even for high energy 7-rays and cosmological 
times, gives an effect that is too small for observable consequences. 

The case (IT3|) is similar to (fHj) . but here we might question if it is acceptable to 
have second order time derivatives on the right hand side of the equation, or we should 
disregard this possibility. 

Similar considerations regarding orders of magnitude and observability hold for the 
parity violating cases (JTSJ), (fTB|). (fT7|). and (fT%|). The last case, with % 4 = — Xs, was 
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obtained in j2|. We notice however that the available observational data indicates that 
Xa — 10 _4 jn], much smaller than the expected value of order one. 

It is worth mentioning that the same observational data can be used to test the 
validity of a recent result that also predicts a rotation of the polarization vector [I]. 
Denoting by D the cosmological distance to the source of the incoming radiation, by 
9gp — X^pk 2 D and 9st — i 2 pL 1 ~ 2a k 2 D the rotation angle of the polarization direction 
obtained by Gambini-Pullin and Sahlmann- Thiemann respectively, and comparing the 
uncertainty of these angles with the observational data, one obtains the following 
restriction for the constants a and L, 




However, these constants must satisfy the following inequalities < a < 5 and 
L » £pjl], which clearly violate the above equation. This would suggest that the 
kinematical states used by the authors to derive their equations cannot be considered 
as semiclassical states. 



3. The interaction between the gravitational and electromagnetic fields 

In this section we first review the standard classical Lagrangian and Hamiltonian 
formulation of coupled Maxwell and gravitational fields. This review is then used to 
argue that, even when we promote the fields to quantum operators, there will always 
be relations that remain valid for the quantum theory. 

The maxwell field on a curved spacetime is given by an exact 2-form 

F = dA, 

with A the Maxwell connection 1-form. Given a local coordinate system with timelike 
and spacelike coordinate vectors e a Q and e™ respectively, we define the electric and 
magnetic fields as = F io and Bi = e^^Fj^, where is the Levi-Civita symbol. 
With these definitions (and assuming that A Q = since we are only concerned with 
propagating waves) the relations between the fields and the potential can be written as 

Ei = -d t A h (19) 

and 

Bi = e ijk djA k . (20) 

Note that these relations do not depend on any metric and are just the coordinate 
components of the above equation. 

Assume now that we give any Lagrangian density C = C(F,g) describing the 
coupling between the Maxwell and gravitational field. By definition, the Maxwell field 
will satisfy 

alF = 0, 
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since it is not part of the Euler-Lagrange equations of motion. The coupling between 
the metric and electromagnetic fields will of course depend on the detailed form of the 
Lagrangian density. In particular, for 

C l = --FA*F (21) 
we obtain the standard field equations 

d*F = 0, 



but we will allow for more general types of couplings that might arise from the 
semiclassical aproximation when we promote the classical fields to quantum operators. 
To construct a Hamiltonian formulation we first define the canonical momentum 
conjugated to A4 



i dC 

71 



d(d t Af) 

and then invert this relation to obtain 

d t Ai = Gi(iT j ,A k ). (22) 

The Legendre transformation 7i = n l Gi — £ then gives the desired Hamiltonian density. 
Note that, by construction, the hamilton equation of motion for will be (j22j) . Note 
also that, regardless of the form of Gj, must satisfy equation (JTHJ). Thus, we 
expect that for a non-standard interaction Hamiltonian, —Ei will not be the conjugate 
momentum to Aj. 

If we now promote the classical fields to quantum operators (A, F, g) we might 
expect several changes in the quantum hamiltonian but the basic definitions should 
remain valid. For example, 

F = dA, 

says that F is still the curvature of A. Thus, the equation 

dF = 0, 

should hold as an identity for A. Moreover, the above relation does not depend on g, 
therefore even if we must regularize the metric operator it should remain unchanged. 
Even if we take expectation values of this relation and assume the state is a direct 
product of a coherent state for the Maxwell field and a "semiclassical state" for the 
gravitational field one should reobtain 



1 class dA c i assy dF c i ass 



at least at a semiclassical approximation. Thus, the plan is as follows: 
(i) We start by impossing F = dA. 
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(ii) We then take the regularized energy density of the electromagnetic field written in 
terms of E± and Bi and write a modified relation between Ei and ir l so that, via 
the hamilton equations, the relation 

Ei = -dtAi 

holds order by order in a perturbation expansion. 

(iii) With the desired relation Ei = —Gi we obtain the second equation of motion for 
7r* and then rewrite the final result as a modified set of equations for the Maxwell 
fields. 



3.1. The phenomenological Hamiltonian and Lagrangian densities 

Following the Gambini-Pullin approach[2|, the phenomenological interaction Hamilto- 
nian density (obtained by taking expectation values of the regularized quantum hamilto- 
nian with coherent states for the Maxwell field and "weave states" for the gravitational 
field) is given by, 

Use = 2 (p 2 + B 2 ) + X l P (E ■ V x E + B ■ V x B) . (23) 

where x is a phenomenological coupling constant that destroys parity invariance. 
The relationship between the fields and the conjugate variables (A, jf) is, by 
assumption, given by 

B = V x A (24) 
E= -vr + xW], (25) 

where we have adopted the vectorial notation for ease of writing and where F[n] is 
a function of n and its derivatives. If we now insert these relations in the above 
Hamiltonian and impose that E = —d t A (which must be true by definition), the 
functional F is determined in a unique way via the Hamilton equation for A, leading to 

E = —7? + 2 X l P V x 7f. (26) 

This allows us to obtain the second equation of motion for jr. Hamilton equations for 
the conjugate variables then read 

d t A = 7T- 2 X l P V x vf , (27) 

d t jf = - V x V x (A + 2 X 1 P V x A). (28) 

Combining these equations and assuming the Coulomb gauge we obtain 

(d 2 -V 2 )A = 0(( X lp) 2 ). (29) 



Alternatively, we can derive Maxwell's equations for the electric and magnetic fields 
yielding 

VxB-d t E = (( X l P ) 2 ) , (30) 
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VxE + d t B = 0(( X l p ) 2 ), (31) 

V • E =0 (( X l p ) 2 ) , (32) 

V ■ B =0. (33) 
Remarks: 



• Note that our field equations are Lorentz invariant and have a standard dispersion 
relation. 

• Starting from the same eq. (}2*3*j) . Gambini and Pullin|2j assumed a different relation 



between the electric field and the momentum, namely E = — k. This assumption 
leads to the following equations of motion for the conjugate pair, 

d t A = - (E + 2 x l p V x E) , (34) 

d t E = Vx(B + 2 x l p SJ x B) . (35) 

In terms of the fields, the equations can be rewritten as 

V xB-d t E = 2 x l p V 2 B, (36) 

V x E + d t B = - 2 x l p V 2 E, (37) 



which, unlike the ones we obtained, are not Lorentz invariant. It is important to 
notice that this breakdown of invariance follows from the assumption that —E and 
A are canonical variables. Although this is the case for the classical theory, it might 
not remain valid when we promote the fields to quantum operators. On the other 
hand, the quantum version of E = —d t A should remain valid since this is just a 
consequence of F — dA. We therefore conjecture that the homogeneous Maxwell 
equations, contained in F = dA, should remain unchanged for the full theory. 

• It is worth mentioning that using a completely different quantization procedure Ellis 
et. al. obtain a generalized set of Maxwell equations that keep the homogeneous 
part unchanged [5]. 

It is also possible to derive, up to linear order in X , the Lagrangian density associated 
to the Hamiltonian ([23)1. via the Legendre transformation £ = it l d t Ai — TC. If we invert 
flUJ) we get 

7? = —E — 2 X l p V x E, (38) 
which leads to the following Lagrangian density 

C EB = ~ (E 2 - B 2 ) + X 1 P (E ■ V x E - B ■ V x B) . (39) 

The Euler-Lagrange equation of motion for the potential can be written in terms of the 
fields as: 

VxB-d t E = -2 X 1 P V x(VxB-d t E)+0 (( X l p ) 2 ) , (40) 

Note also that in principle the above equations could give more solutions than plane 
waves since they are second order PDE's. Thus, there seems to be more solutions to the 
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above equations than Hamiltons equations (|3U|) . since the ones obtained directly from 
the Hamiltonian do not contain terms linear in \. The extra terms appearing in (J4~U|) 
arise when we invert the relation between E and tt. We can see that both equations 
contain the same solutions, in the linear approximation, by the following consideration. 
Let C = V x B — dfE. The solution to Hamilton equations correspond to C — 0. We 
now search for non trivial solutions to 

2 X 1 P V xC + C = 0, (41) 

Note that if C satisfies (14 lj) . then it is also a solution of 

C+(2 X 1 P ) 2 V 2 C = 0, (42) 

and this implies that solutions of ()41|) have a "monochromatic" (spatial) Fourier 
spectrum with k = (2%/p) -1 . Such contributions with wavelength of the order of the 
Planck length must be absent if the low energy approximation (characteristic lengths 
much larger than l p ) is to be consistent. 

Assuming this cut off is implemented, the only solution consistent with our 
approximation is V x B - d t E = 0. 

Using a somewhat different approach, L. UrrutiajB] developed a Lagrangian 
formulation that leads to equations (jHBl ET7|) . Starting from (j2~Tj) and defining E by 
Eq. (|H4*|) and B as the curl of A, Urrutia obtains a Lagrangian density for A such that 
the corresponding Euler-Lagrange equations yield (JHEl IHTj) . This Lagrangian, however, 
is a non-local functional of A. The fact that a Lorentz violating dispersion relation 
uj± = ^Jk 2 T- ^xl P k 3 — 1^1(1 =F 2x^ P |fc|) is found for the set of equations (|HH| ITTj) is 
consistent with a recent result advocating that non-local Lagrangians can generate non- 
Lorentz-invariant dispersion lawsjHj- The non-local Lagrangian arose, however, from an 
assumed non-local relation between E and A which will not be true if one follows our 
assumption. 

Unlike the approach of UrrutiajB] the Lagrangian density that yields eqs. tjHOJHlj) is a 
local functional of the potential. This result is in agreement with recent work of J. Bros 
and H. Epstein ^0] which states that microcausality and energy positivity in all frames 
imply Lorentz invariance of dispersion laws. (Within the domain of our assumption 
the energy spectrum will be positive definite since the Hamiltonian can be written as 
H S b = \{e + xhN * E) + | (B + x/pV x b) plus terms of order x 2 )- 

It would appear that the issue of whether or not Lorentz invariance is broken 
depends on the specific relation between the fields and the potential. The local 
Lagrangian constructed from the potential and its derivatives yields Lorentz invariant 
dispersion laws, whereas the non-local Lagrangian breaks the invariance. Following our 
assumption that F = dA we conjecture that the quantum fields will be local functionals 
of A. 

In principle, if the Hamiltonian were given as a perturbation expansion in x, it 
would be possible to introduce higher order corrections to the relationship between the 
electric field and conjugate momentum, eq. (J2H)) . so that Faraday's law is preserved at 
each order in that expansion. 
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4. Final comments and conclusions 



We first summarize the main results of this work. In Section 2 we concentrate on the 
phenomenological field equations for Maxwell fields that could arise from the interaction 
with the gravitational field in the semiclassical approximation. We show that, up 
to linear order, there is just one acceptable way to obtain a non-standard dispersion 
relation, namely, field equations with parity violating coupling constants[2.. 

Is the universe living in a state of definite parity? Voting for the affirmative, the 
cosmological nature and detailed temporal structure of gamma ray bursts could serve as 
a tool to observe the predicted deviations from the standard dispersion relations. As was 
shown on a previous workJS] the phenomenological constant is rather small and future 
observations could be useful to decide if indeed it is non- vanishing. However, using the 
available data we can essentially rule out the Sahlmann-Thiemann construction. 

In view of the present observational data for cosmological sources and the 
ambiguities of the available models for the interaction of Maxwell fields with a 
semiclassical space-time it is fair to ask ourselves whether or not we should still have 
faith in Lorentz invariance. In Section 3 we show that a very natural assumption leads 
to Lorentz invariant field equations. Our conclusion is that we have no reason to believe 
that a quantum theory of gravity would change the invariance and thus conjecture that 
Lorentz invariance will still hold in the quantum interaction between gravitons and 
photons. 
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